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DURATION: 

ACCURACY: 
DESCRIPTION: 



UNIVERSITY OF ILLINOIS 

DIGITAL COMPUTER 

AUXILIARY 

LIBRARY ROUTIN E E 9 - 308 

Closed (Nevton -Codes) and Open (Steffensen) Quadrature 

(Numerical Integration for tabulated function values) 

Closed subroutine 

15 ■+ n (Closed Quadrature) 

13 + n (Open Quadrature) 

(n defined below) 



2.9 + l.^n milliseconds (Closed Quadrature) 

2.9 + lA(n-2) milliseconds (Open Quadrature) 

+ 2 + truncation error (see below) 

This routine calculates an approximation to the integral, 

using tabulated function values, 



a 



to Z $ Z N k f k' 

D k=0 * * 



where 



f k - f(a + k %ai). 



n 



USE: 



The closed and open quadrature cases are distinguished by 
the fact that N Q and N are zero for open quadrature, and 
non-zero for closed quadrature. This means that an approxima- 
tion to the integral can be computed when the terminal function 
values are not known. The number n is the number of sub- 
intervals into which the original interval (b-a) is divided; 
thus there are n + 1 function values in the closed case, and 
n - 1 function values In the open case. The numbers D and N, 
are chosen to give optimum results for equally spaced values of 
x within the region of integration. 

To use this routine, the programmer will determine the number 
of intervals in his region of integration, and the type of 
quadrature desired. The coefficients appropriate to his choice 
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are then copied from the second part of the library tape 
immediately following the program. The routine is entered 
with the orders 



TRUNCATION ERROR: 



p 


— tF 
50 pF 


p+1 


26 xF 



where x is the location of the first word of this routine. 

The program parameter t is the location of the first of the 

tabulated function values f ,; ■ in the closed case, the function 

value f, will be in location t + k. In the open case, f. 

is to be at t, f ' at t + 1, ..., and f^ at t + n-2. (Because 

N = N = in the onen case, the function values f n and f 
On' u n 

are neither needed nor used by the routine)* 

When control is returned to the right-hand order in p + 1, 

the computed integral will be in both the accumulator and 

quotient registers. 

The maximum error introduced in the integration by using only 

a finite number of function values is given approximately 

by the following table: 



n=number 

of intervals 



1 
2 

3 
k 

5 

6 

7 
8 

9 

10 
11 
12 



Closed Case 



(z) 



-3.3x10 f^ ' (z) 
-l.bxlO f K ' (z) 
-5.2xlO _T f (6) (z) 
-3.0XI0- 7 f (6) (z) 
-6.4X10- 10 f (8) (z) 
-4.0X10" 10 f (8) (z) 
-5.9X10- 13 f (l0) (z) 
-3.8X10- 15 f (l0) (z) 
A.2X10' 16 f (l2) (z) 
-2.7X10- 16 f (l2) (,)■ 
-2.2X10" 19 f (l4) (z) 



Open Case 



T.W* f V (Z) 
2.2X10" 4 f (4) (z) 
l.lxlO" 6 f (6) (z) 
7.1txl0" 7 f (6) (z) 
2.^10" 9 f (8) (z) 
1.5X10" 9 f (8) (z) 
2.TX10- 12 f (l0) (z) 
2.0X10" 12 f (l0) (z) 
2.5xlO" 15 f (l2) (z) 
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NOTE 1: 



NOTE 2: 



REFERENCES: 



Here, f ^ (z) is the pth derivative of the function being 

integrated, evaluated at some point z in the interval (a,b). 

This routine is independent of a_ linear change of variable, 

so that the limits of integration a and b need not be specified, 

If the programmer does not need to use function values of 

equally spaced points (for example, the function values are 

evaluated by an auxiliary subroutine), then greater accuracy 

may be obtained by using Gaussian Quadrature (see library 

routines pertaining to Gaussian Quadrature). 

B. P. Moors, "Valeur Approximative d'une Integrale Defini" 

J, F. Steffensen, "Interpolation" 
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OOK 









K5 F 








42 11L 


Plant link 




1 


46 5L 








1+1 F 


Plant function value address, clear sum box 




2 


L5 3L 








42 4L 


Set coefficient address 




3 


23 4L 








00 l4L 


enter loop 




4 


S5 F 
50 F 






5 


7^ F 
L4 f 






6 


40 F 










Step function value list address 




7 


L4 6l 
46 5L 






8 


F5 4l 








42 4l 


Step coefficient address 




9 


LO 12L 








36 4L 


test for end 




10 


L5 F 
66 13L 






11 


S5 F 
22 F 


5 Zf k N k 




12 


35 F 








50 L 


End constant 




13 


00 F 








00 F 


D 




14 


00 F 








00 F 


Coefficients N Q (closed), N, (open) 




9 




• • • 

to be 




• 




copied 




• 




from 
library 






00 F 


tape. 






00 F 


N n N „-2 



